3:17 Integral Transforms and Special Functions gegpoly 

Integral Transforms and Special Functions 
Vol. 00, No. 00, January 2009, 1-12 



RESEARCH ARTICLE 

On a generalization of the generating function for Gegenbauer 

polynomials 

Howard S. CohP* 

'^Applied and Computational Mathematics Division, National Institute of Standards and 
Technology, Gaithersburg, Maryland, U.S.A. 

{v3.5 released October 2008) 

A generalization of the generating function for Gegenbauer polynomials is introduced whose 
coefficients are given in terms of associated Legendre functions of the second kind. We discuss 
how our expansion represents a generalization of several previously derived formulae such as 
Heine's formula and Heine's reciprocal square-root identity. We also show how this expan- 
sion can be used to compute hyperspherical harmonic expansions for power-law fundamental 
solutions of the polyharmonic equation. 
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1. Introduction 

Gegenbauer polynomials C^{x) are given as the coefficients of p" for the generat- 
ing function (1 + /o^ — 2px)~'^ . The study of these polynomials was pioneered in a 
series of papers by Leopold Gegenbauer (Gegenbauer (1874,1877,1884,1888,1893) 
(l3l-[l7|). The main result which this paper relies upon is THEOREM l2.1l below. This 
theorem gives a generalized expansion over Gegenbauer polynomials Cn{x) of the 
algebraic function z i— t- {z — x)~" . Our proof is combinatoric in nature and has 
great potential for proving new expansion formulae which generalize generating 
functions. Our methodology can in principle be applied to any generating function 
for hypergeometric orthogonal polynomials, of which there are many (see for in- 
stance Srivastava & Manocha (1984) Erdelyi et al. (1981) [lH). The concept 
of the proof is to start with a generating function and use a connection formula to 
express the orthogonal polynomial as a finite series in polynomials of the same type 
with different parameters. The resulting formulae will then produce new expansions 
for the polynomials which result from a limiting process, e.g., Legendre polynomi- 
als and Chebyshev polynomials of the first and second kind. Connection formulae 
for classical orthogonal polynomials and their (/-extensions are well-known (see Is- 
mail (2005) [20]). In this paper we applied this method of proof to the generating 
function for Gegenbauer polynomials. 

This paper is organized as follows. In Section[2]we derive a complex generalization 
of the generating function for Gegenbauer polynomials. In Section[3]we discuss how 
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2 Howard S. Cohl 

our complex generalization reduces to previously derived expressions and leads to 
extensions in appropriate limits. In Section |4] we use our complex expansion to 
generalize a formula originally developed by Sack (1964) (23 ] on R^, to compute an 
expansion in terms of Gegenbauer polynomials for complex powers of the distance 
between two points on a d-dimensional Euclidean space for d>2. 

Throughout this paper we rely on the following definitions. For ai, a2, as, . . . G C, 
if i,j G Z and j <i then X]n=i ~ ^ ^^"^ nn=i'^" ~ where C represents the 
complex numbers. The set of natural numbers is given by N := {1, 2, 3, . . .}, the 
set No := {0,1,2,...} = N U {0}, and the set Z := {0, ±1, ±2, . . .}. The sets Q 
and R represent the rational and real numbers respectively. 



2. Generalization of the generating function for Gegenbauer polynomials 

We present the following generalization of the generating function for Gegenbauer 
polynomials whose coefficients are given in terms of associated Legendre functions 
of the second kind. 

Theorem 2.1 LetueC\ -Nq, G (-1/2, oo) \ {0} and z e C \ (-co, 1] on any 
ellipse with foci at ±1 with x in the interior of that ellipse. Then 



1 2^+V2r('M)e*'^(''-'"+i/2) ^ 



If one substitutes z = (1 + p^)/(2p) in ([T]) with < |p| < 1, then one obtains an 
alternate expression with x G [—1, 1], 



r(^)e 



iTv(p-u+l/2) 



(1 + p2 _ 2pxy V^r(z^)p^+V2(l _ p2y^f^-l/2 

oo / 1 _i_ 2 \ 

><E(-+/^)Q:;ri/2Rf K^n^w- (2) 

n=0 ^ ^ ^ 



One can see that by replacing = /i in ([2|), and using (8.6.11) in Abramowitz &: 
Stegun (1972) [l[, that these formulae are generalizations of the generating function 
for Gegenbauer polynomials (first occurence in Gegenbauer (1874) [13]) 



where /3 G C with 
Olver et al. (2010) 

by 



(l + ,2_,,,). =EC(x).", (3) 

£l < 1 and u G (-1/2, oo) \ {0} (see for instance (18.12.4) in 
23]). The Gegenbauer polynomials : C — t- C can be defined 



C'^.ix) := (-n, n + 2.; . + 1; i^) , (4) 

where n G Nq, G (-1/2, oo) \ {0}, and 2F1 : x (C \ -Nq) x {z G C : 
\z\ < 1} — 7- C, the Gauss hypergeometric function, can be defined in terms of the 
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following infinite series 

2Fi{a,b;c;z) ■.= }^—- — (5) 

n=0 

(see (2.1.5) in Andrews, Askey & Roy 1999), and elsewhere by analytic continua- 
tion. The Pochhammer symbol (rising factorial) {■)n : C — )• C is defined by 



{z)n :=Y[iz + i-l), 

i=l 

where n S Nq. For the Gegenbauer polynomials C^{x), we refer to n and as the 
degree and order respectively. 

Proof Consider the generating function for Gegenbauer polynomials ([3]). The con- 
nection relation which expresses a Gegenbauer polynomial with order v as a sum 
over Gegenbauer polynomials with order is given by 

(2t/)„ ^ {i^ + k + l)n-ki2i^ + n)kifi + ^)kTi2n + k) 
"^""^ " (^ + l)n ^0 (n - k)\ (2^), T{2^, + 2k) 

^^^K l^ + k + l,2^. + 2k + l ^'V^^'(-)- (6) 

This connection relation can be derived by starting with Theorem 9.1.1 in Ismail 
(2005) [2Q] combined with (see for instance (18.7.1) in Olver et al. (2010) [13]) 

c(x) = -^pr^/^-^-^/^)(x), (7) 

[1'+ i)n 



i.e., the Gegenbauer polynomials are given as symmetric Jacobi polynomials. The 
Jacobi polynomials Pj"'^^ : C — )• C can be defined as (see for instance (18.5.7) in 
Olver et al. (2010) [23]) 

Pf^ (x) := ^^^p^ 2F1 (^-n, n + a + (3 + l;a + l; , (8) 



where n G Nq, and a,p > -1 (see Table 18.3.1 in Olver et al. (2010) [23]). The 
generalized hypergeometric function 3F2 : x (C\— Nq)^ x {z G C : |2;| < 1} — )• C 
can be defined in terms of the following infinite series 

/ai, 02,03 \ {ai)n{a2)n{a3)n z"' 

b..bo ^=Z^ 



61,62 ' y ■ (&l)"(&2)„ n! 

n=0 



If we replace the Gegenbauer polynomial in the generating function ([3]) using the 
connection relation ([6]), we obtain a double summation expression over k and n. 
By reversing the order of the summations (justification by Tannery's theorem) and 
shifting the n-index by k, we obtain after making some reductions and simplifica- 
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tions, the following double-summation representation 



k 



{l + p'^-2pxY 22'^-ir(z.)^ '2'^kT{u + k + \)T{p + k + l) 

^ T{2v + 2k) /-n,n + 2k + 2iy,fi + k+l \ , . 

^ n\ ^ + k + ^,2fi + 2k + l 'V' 



X 

n=0 



The 3-F2 generalized hypergeometric function appearing the above formula may be 
simplified using Watson's sum 



a, b, c 
Ka + 6+ l),2c' 



3F2 1. . : ' . ;1 



r{u + k + l)r{fi + k + i) 



_ V^T (c + ^) r (g + 6 + 1)) r (c + ^ (1 - g - b)) 

"r(i(g + i))r(i(6 + i))r(c + i(i-a))r(c + i(i-6))' 

where Re(2c — a — b) > —1 (see for instance (16.4.6) in Olver et al. [23]), therefore 

f -n,n + 2k + 2v,p + k + \ \ 
^ ^\ y + k+\,2^i + 2k + l ' ) 

^ ^T{ji-v + l) 

r (if^) T{v + k + ^)T{p + k + l + ^)T{p-u + l-^) 

for Re(/U — v) > —1. By inserting (llOp in ([9]), it follows that 

1 _ ^r(^)r(/.-. + i) ^ 

(1 + p2 _ 2pxY ~ 22-ir(z.) ^ >^ ^ > 22^= 



, (10) 



oo 

X 



fe=0 

p"r(2z/ + 2k + n) 



^o"'r(i^)r(z. + fc + n±i)r(/. + A: + i + f)r(/.-^ + i-f] 

It is straightforward to show using ([5]) and 

^n r(z) 



r(z-n) = (-l)^ 



(-^ + l)n' 



for n G No and z G C \ — Nq, and the duplication formula (i.e., (5.5.5) in Olver et 
al. (2010) [13) 



222-1 / -y- 

r(2.) = -^r(.)r . + -), 



January 18, 2013 3:17 Integral Transforms and Special Functions gegpoly 



Integral Transforms and Special Functions 

provided 2z — Nq, that 



E 

n=0 



/)"r(2i/ + 2k + n) 



n\T (i^) r (zv + A: + ^) r (/i + A: + 1 + §) r (/i - + 1 - f ) 



so therefore 



:-n^i±^2Fl{^ + t,^-K^ + *: + l;p2). (11) 



Finahy utiHzing the quadratic transformation of the hypergeometric function 



' a a -\- \ 4lZ 
2Fi{a,b;a-b+l;z) = (1 + z)-'^2i^i ( ;a-6+l; 



2' 2 ' '(z + 1 



\2 



for |z| < 1 (see (3.1.9) in Andrews, Askey & Roy (1999) [1]), combined with the 
definition of the associated Legendre function of the second kind : C\(— oo, 1] — )• 
C in terms of the Gauss hypergeometric function 



g/.(^) V^e^^'^Tiu + ^ + l)(z2 _ 



,V + ^ + 2 u + ii + l 3 1, 



for |z| > 1 and 1/ + /i + 1 ^ -Nq (cf. Section 14.21 and (14.3.7) in Olver et al. 
(2010) [ii]), one can show that 



2F1 [u + k, u - fi; fi + k + 1; p- 



2^ 



r (^ + A; + 1) e^-(/^--+i/2) fl + p' 



v^r(i/ + k)p^'+^+y^{l - p2)z.-M-i/2 ^fc+M-1/2 \^ 2p 



which when used in ()11|) produces ([2]). Since the Gegenbauer polynomial is just a 
symmetric Jacobi polynomial ([7]), through THEOREM 9.1.1 in Szego (1959) (26| 
(Expansion of an analytic function in a Jacobi series) , since : C — )• C defined by 
fz{x) := {z — x)~'^ is analytic in [—1, 1], then the above expansion in Gegenbauer 
polynomials is convergent if the point z £ C lies on any ellipse with foci at ±1 
and X can lie on any point interior to that ellipse. H 
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3. Generalizations, Extensions and Applications 

By considering in ([2]) the substitution v = d/2 — 1 and the map v i— t- — one 
obtains the formula 



P 



oo / 1 I 2 \ 

X j:(2n + d - 2)Q:^; ( iif ) C^J^-Hx). 

n=0 ^ ^ 

This formula generalizes (9.9.2) in Andrews, Askey & Roy (1999) f^. 

By taking the limit as /i — )• 1/2 in ([T|), one obtains a general result is an expansion 
over Legendre polynomials, namely 

= ^-YUa^ E(2- + l)Qr^(^)^'n(x), (13) 

using 

Pn{x) = Cy2(^), (14) 

which is clear by comparing (cf. (18.7.9) of Olver et al. (2010) and (jj]) or ([8])) 

Pn{x):=2Fi(-n,n + l-l;^^\, (15) 



and (H]). If one takes = 1 in (jl3p then one has an expansion of the Cauch\ 
denominator which generalizes Heine's formula (see for instance Olver (1997) [22: 
Ex. 13.1]; Heine (1878) 0, p. 78]) 



^ oo 

— = Y,{'^n + l)Qn{z)Pn{x). 



z — X 

ra=0 

By taking the limit as /i — )• 1 in ([T|), one obtains a general result which is an 
expansion over Chebyshev polynomials of the second kind, namely 

(7^ = 0Fr(.)(.2- 1)^/2-3/4 + l)Q„+l/2(-)f^n(x), (16) 

V \ / \ / n=0 

using (18.7.4) in Olver et al. (2010) [1^, Unix) = C^(x). If one considers the case 
= 1 in (jl6[) then the associated Legendre function of the second kind reduces 
to an elementary function through (8.6.11) in Abramowitz & Stegun (1972) [l[, 
namely 



Unix) 



+1 



By taking the limit as — t- 1 in ([T]), one produces the Gegenbauer ex pan sion of 
the Cauchy denominator given in Durand, Fishbane & Simmons (1976) [lOj, (7.2)], 
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namely 

1 2^+1/2 _ 



Z — X 

^ n=0 

Using (2.4) therein, the associated Legendre function of the second kind is converted 
to the Gegenbauer function of the second kind. 
If one considers the hmit as /i — )• in ([1]) using 

lim ^C/^(x) = 6„r„(x) 

(see for instance (6.4.13) in Andrews, Askey & Roy (1999) 0]), where T„ : C — )• C is 
the Chebyshev polynomial of the first kind defined as (see Section 5.7.2 in Magnus, 
Oberhettinger & Soni (1966) [IH) 



2' 2 

and €n = 2 — 6nfi is the Neumann factor, commonly appearing in Fourier cosine 
series, then one obtains 

/T„-i7r(!/-l/2)/^2 i\-i//2+l/4 °° 



(z - x)'^ V vr r(zy) 

^ ^ ^ ^ n=0 

The result (I17p . which is a generalization of Heine's reciprocal square-root identity 
(see Heine (1881) 0, p. 286]; Cohl & Tohline (1999) @, (A5)]). Polynomials in 
{z — x) also naturally arise by consideringthe limit ly ^ n € — Nq. This limit is 
given in (4.4) of Cohl & Dominici (2010) [71], namely 

(. - xr = ^{-lr'Jliz' - ly'^^''' t enT.i-)j^Qt\%iz), (18) 

n=o 

for q G Nq. 

Note that all of the above formulae are restricted by the convergence criterion given 
by Theorem 9.1.1 in Szego (1959) [26^ (Expansion of an analytic function in a 
Jacobi series), i.e., since the functions on the left-hand side are analytic in [—1,1], 
then the expansion formulae are convergent if the point z £ C lies on any ellipse 
with foci at ±1 then x can lie on any point interior to that ellipse. Except of course 
(jlSp which converges for all points z,x G C since the function on the left-hand side 
is entire. 

An interesting extension of the results presented in this paper, originally uploaded 
to arXiv in Cohl (2011) [6] have been obtained recently in Szmytkowski (2011) [2^ . 
to obtain formulas such as 



A. ^ ^„+^-i/2W^nW 2M-i/2r(/x + i)r (1 - Z.) 



if-l<x<t<l, 
{x - t)-"-^/'^ if - 1< t < X < 1, 
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and 



n=0 



{t - x)-"-^/^ if - 1< X < t < 1, 

ty-^/^ cos[n{u + i)] if - 1< t < x < 1, 



[X 



where Re^ > —1/2, Rei/ < 1/2 and Pu,Q,u : (—1,1) — )■ C are Ferrers functions 
(associated Legendre functions on-the-cut) of the first and second kind. The Ferrers 
functions of the first and second kind can be defined using Olver et al. (2010) 0, 
(14.3.11-12)]. 



4. Expansion of a power-law fundamental solution of the polyharmonic 
equation 



A fundamental solution for the polyharmonic equation on Euclidean space R is a 



function Qi : (R*^ x R"^) \ {(x, x) : x G R'^} — )• R which satisfies the equation 



(-A)'=gf(x,x') = 5(x-x'), (19) 
where A : CP(R'^) — ^ CP^^(R'^), p > 2, is the Laplacian operator on R*^ defined by 



d q2 
j=l ' 

X = (xi, . . . ,Xrf),x' = {x[, . . . ,x'^) G R'', and 6 is the Dirac delta function. Note 
that we introduce a minus sign into the equations where the Laplacian is used, such 
as in p9p. to make the resulting operator positive. By Euclidean space R*^, we mean 
the normed vector space given by the pair (R^, II • II), where || • || : R'^ ^ [0, oo) is 

the Euclidean norm on R'^ defined by ||x|| := \J x\ + ■ ■ ■ + x^, with inner product 

(•, •) : R"' X R"' ^ R defined as 

d 

(x,x') := J^x.x^ (20) 

i=l 

Then R*^ is a C°° Riemannian manifold with Riemannian metric induced from 
the inner product ([20]). Set S^"^ = {x G R'^ : (x,x) = 1}, then S'^-\ the {d - 1)- 
dimensional unit hypersphere, is a regular submanifold of R^ and a Riemannian 
manifold with Riemannian metric induced from that on R"^. 



Theorem 4.1 Let d,k eN. Define 

^_l-jk+d/2+l ||x — x'lp'^"'^ 



gf(x,x') : = 



(log ||x - x'll - /3p 



{k-iy. {k- d/2)\ 22fc-%'=f/2 ^^"^ "^^ " ^P^'^' 

if d even, k > d/2, 

r(d/2- A:)||x-x'|P*^-'^ 

— -- — j7- otherwise, 

(/c - 1)! 22fc7r°'/2 
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where p = k - d/2, (3p^d ^ Q is defined as j3p^d ■= \ [Hp + Hd/2+p-i - -ffd/2-1] 
with Hj £ Q being the jth harmonic number 



^ 1 



I 

i=l 



then Q'^ is a fundamental solution for {—A)^ on Euclidean space R'^. 
Proof See Cohl (2010) |4.] and Boyling (1996) [1]. 

Consider the following functions g^, : (R'^ x R'^) \ {(x, x) : x € R'^} — )• R 
defined for d odd and for d even with k < d/2 — 1 as a power-law, namely 

Sf(x,x'):=||x-x'f^-^ (21) 

and for d even, k > d/2, with logarithmic behavior as 

[f(x,x') := ||x-x'fP(log||x-x'||-/3p,,), 

with p = k — d/2. By Theorem 14. II we see that the functions gf. and if equal real 
non-zero constant multiples of for appropriate parameters. Therefore by (|19|) . 

and if are fundamental solutions of the polyharmonic equation for appropriate 
parameters. In this paper, we only consider functions with power-law behavior, 
although in future publications we will consider the logarithmic case (see Cohl 
(2012) [| for the relevant Fourier expansions). 

Now we consider the set of hyper spherical coordinate systems which parametrize 
points on R'^. The Euclidean distance between two points represented in these 
coordinate systems is given by 

||x — x'll = V2rr' [z — cos 7]^''^ , 

where the toroidal parameter z € (l,oo), (2.6) in Cohl et al. (2001) is given by 



z :- 



j,2 _j_ ^'2 

2rr' ' 



and the separation angle 7 G [0, vr] is given through 



(x,x') , , 

cos7 = ^^, (22) 



where r, r' G [0, 00) are defined such that r := ||x|| and r' := ||x'||. We will use these 
quantities to derive Gegenbauer expansions of power-law fundamental solutions for 
powers of the Laplacian g^ (i2T]) represented in hyperspherical coordinates. 

Corollary 4.2 For d e {3,4,5, . . .}, u £ C, x, x' G R'' with r = ||x||, r' = ||x'||, 
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and COS 7 = (x, x')/(rr'), the following formula holds 

°° , , , /^2 I ^/2^ 



ix-x'ir 



(2n + d - 2) QiV(r?//2 ( SJ;^ 7), (23) 



X 

ra=0 



w/iere := "^^^{r^r']. 

Note that (|23p is seen to be a generalization of Laplace's expansion on (see for 
instance Sack (1964) [H]) 



1 v—^ r" 



X — X' 

' " n=0 > 



which is demonstrated by utilizing (|14p and simplifying the associated Legendre 

'2 

1/2 



function of the second kind in (f23]l through Q\i2 ■ C \ (— oo, 1] — )• C defined such 



that 

Q%iz) = ^^{z^ - 1)-V4 

(cf. (8.6.10) in Abramowitz & Stegun (1972) fl*] and (fT2D). 

The addition theorem for hyperspherical harmonics, which generalizes 

P„(C0S7) = -^ V Yn,m{Sl)Yn,m{^'), (24) 

2n + 1 ^-^ 



m=—n 



where Pnix) is the Legendre polynomial of degree n £ Nq, for d = 3, is given by 

Y^^i^)W^) = 2^1/1 (f! 2) + " 2)^n/'-'(cos 7), (25) 

where K stands for a set of {d — 2)-quantum numbers identifying degenerate har- 
monics for a given value of n and d, and 7 is the separation angle (j22p . The functions 
1^ : S'^~^ — )• C are the normalized hyperspherical harmonics, and Yn,k : — )• C 
are the normalized spherical harmonics for d = 3. Note that x, x' € S'^"^, are the 
vectors of unit length in the direction of x, x' G R*^ respectively. For a proof of the 
addition theorem for hyperspherical harmonics ()25p . see Wen &: Avery (1985) [sO] 
and for a relevant discussion, see Section 10.2.1 in Fano & Rau (1996) [12]. The 
correspondence between (p^ and (j25|) arises from dH) and (fT5]) namely p^ . 

One can use the addition theorem for hyperspherical harmonics to expand a fun- 
damental solution of the polyharmonic equation on R*^. Through the use of the 
addition theorem for hyperspherical harmonics we see that the Gegenbauer poly- 
nomials Cn^'^ ^(cos7) are hyperspherical harmonics when regarded as a function 
of X only (see Vilenkin (1968) [1^). Normalization of the hyperspherical harmonics 
is accomplished through the following integral 



/ (X)y/(X)di^ = 1, 



January 18, 2013 3:17 Integral Transforms and Special Functions gegpoly 



Integral Transforms and Special Functions 11 

where dO, is the Riemannian volume measure on S"^"^. The degeneracy, i.e., number 
of Unearly independent solutions for a particular value of n and d, for the space of 
hyper spherical harmonics is given by 

Ad — 3 + n)\ 



(see (9.2.11) in Vilenkin (1968) [29l|). The total number of linearly independent 
solutions (j26p can be determined by counting the total number of terms in the 
sum over K in (j25p . Note that this formula (j26p reduces to the standard result in 
d = 3 with a degeneracy given by 2n + 1 and in d = 4 with a degeneracy given by 
(n + l)2. 

One can show the consistency of Corollary 14.21 with the result for d = 2 given 
by 



1^ _ ^ ^_ ■ v>-r|)(^+^/^ ^ ,^(^_^,) (,+i)/2 / + 

' " r(-r./2) ^^-1/2 1 2rr' 



where v G C \ {0,2,4, . . .}, by considering the limit as /i — >• in ([T]) (see pT|) 
above). These expansions are useful in that they allow one to perform azimuthal 
Fourier and Gegenbauer polynomial analysis for power-law fundamental solutions 
of the polyharmonic equation on R'^. 



5. Conclusion 

In this paper, we introduced a generalization of the generating function for Gegen- 
bauer polynomials which allows one to expand arbitrary powers of the distance 
between two points on d-dimensional Euclidean space R"^ in terms of hyperspher- 
ical harmonics. This result has already found physical applications such as in 
Szmytkowski (2011) [13] who uses this result to obtain a solution of the momentum- 
space Schrodinger equation for bound states of the d-dimensional Coulomb prob- 
lem. The Gegenbauer expansions presented in this paper can be used in conjunction 
with corresponding Fourier expansions (Cohl & Dominici (2010) 0]) to generate 
infinite sequences of addition theorems for the Fourier coefficients (see Cohl (2010) 
Q]) of these expansions. In future publications, we will present some of these addi- 
tion theorems as well as extensions related to Fourier and Gegenbauer expansions 
for logarithmic fundamental solutions of the polyharmonic equation. 
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